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Abstract 

Exact wavefunctions for N non-Abelian Chern-Simons (NACS) particles are ob- 
tained by the ladder operator approach. The same method has previously been 
applied to construct exact wavefunctions for multi-anyon systems. The two distinct 
base states of the NACS particles that we use are multi-valued and are defined in 
terms of path ordered line integrals. Only strings of operators that preserve the mon- 
odromy properties of these base states are allowed to act on them to generate new 
states. 
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1. Introduction 



It is now well-known that anyons — particles with arbitrary statistics can exist 
in 2+1 dimensions.' 11 Owing to the topological gauge potential, even non-interacting 
two-anyon states are not (symmetrized) tensor products of single anyon states.' 21 Thus 
the quantum mechanics of systems of many anyons presents a challenge to field the- 
orists. In general, the center of mass motion of the system is not sensitive to the 
statistics and can be factored out. For two anyons, the relative coordinates present us 
with a one-body problem which can be solved in many cases. For N anyons, there are 
N— 1 relative coordinates, whereas there are N(N — 1)/2 pairs of particles. These two 
numbers match only when N = 2. For N > 2, the various pair-separation coordinates 
(in terms of which the statistical gauge potential is easy to write down) are not inde- 
pendent of each other. For this reason, not a single three-anyon system has been com- 
pletely solved. The many-anyon system has been studied in the mean field approach.' 41 
Other methods that have been employed include the semi-classical approximation 151 , 
numerical studies' 61 , perturbative analysis from the bosonic or fermionic ends' 71 , and 
the most interesting of all, the ladder operator approach: In the last couple of years, 
a substantial subset of the exact multi-anyon wavefunctions in a magnetic field has 
been found with this systematic analysis.' 81 This soon generalized to free anyons' 91 and 
anyons of multi-species.' 101 While a lot of states are still missing, all the states in the 
lowest Landau level are obtained by this method. 

In this paper, we apply the ladder operator approach to non-Abelian Chern- 
Simons particles which may be regarded as a generalisation of anyons.' 11 ' In the gauge 
that the Hamiltonian is a free Hamiltonian, the wavefunctions (which have more 
than one component) are multivalued with non-trivial monodromy properties given 

[12] 

by a monodromy matrix. By introducing statistical gauge potentials, one has the 
liberty to work with single-valued wavefunctions. However, since we find multivalued 
wavefunctions convenient to work with, we will stick to them in the rest of this paper. 

In section 2, we review the ladder operator formalism as applied to anyons. This 
helps to highlight the differences between the cases of anyons and non-Abelian C-S 
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particles, the object under study in section 3. In particular, of the operators used for 
anyons, only a subclass of operators, which preserve the monodromy properties of the 
wavef unctions, are allowed to act on the C-S particles. Nonetheless, our wavefunctions 
do cover the lowest Landau level. As an application of our formalism, we compute the 
second virial coefficient of NACS particles. The same set of ladder operators apply 
to free NACS particles with minor modifications. We also consider systems of multi- 
species NACS particles. Finally, the relevance of our work to systems of vortices of 
finite gauge groups is also discussed. 



2. Anyons 

The Hamiltonian for N anyons with charge e and mass m moving on a plane with 
a constant magnetic field B (perpendicular to the plane) is given by 



N I 



a=l 



where the external gauge field A 1 = — ^Be^x^ in the symmetric gauge and the sta- 
tistical gauge potential 



J j 



a l a (x 1} ...,x N ) = vJ2^T^ 12- ( 2 ) 

| X Q x /3| 



By a singular gauge transformation, we can remove a Q from the Hamiltonian at 
the expense of using multi-valued wavefunction 

Vw(xi, • • • ,Xjv) = ex p( w ^2 0q/3 Wo/d(xi, • • • ,xjv). (3) 
Using the complex notation z = x 1 + ix 2 ,z = x 1 — ix 2 ,d = -^,0 = the gauge 
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[8] 

transformed Hamiltonian becomes 



N < ■ .2 D 2 

W? ) - ^J, (4) 



/ 2 - 

# = V <9«<9a + 

^ V m 

a=l N 



e 2 £? 2 , , 9 \ e£ 



8m / 2m 



where J is the angular momentum operator in the singular gauge 

N 

J = ^2(z a d a - z a d a ). (5) 

a=l 

Its eigenvalues are shifted from those in the symmetric gauge by a constant \vN(N — 
1). (See (10. a).) It is convenient to extract a factor exp(— \eB X^=i ka| 2 ) from the 
wavefunction. Then the eigenvalue problem becomes 

H*P = {E- \nuj)^, (6.a) 

Jijj = jxj>, (6.6) 
(with u = where the new Hamiltonian i7 and wavefunctions ^ are defined by 

N 

I v _ > / ; _ \ 

(7.a) 



# = ( -— d a d a + —Z a da] , 
' \ m ml 

a=l N ' 



^ = expf |^ a | 2 j^. (7.6) 

Note that the ground state energy is shifted by \Nu. We impose two physical re- 
quirements for the wavefunctions. Firstly, they must vanish at points of coincidences 
if v 7^ due to the centrifugal potentials (hard-core requirement). Secondly, they 
form Abelian representations of the braid group. 
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Now we introduce the operators 



da — Z a ^j^d a , &a — d a , (8.d) 



t 2 - 

b{ x = z a - ^g d a, b a = d a , (8.6) 

which satisfy [a a , at] = [b a , ot] = <5 a/ 3, all other commutators being zero. With respect 
to these operators, the Hamiltonian H in (7. a) and the angular momentum J in (7.b) 
can be rewritten as 

N t 

a=l 



H = u^2ata a , (9. a) 



N t t 

J = ^2(blb a - a\a a ). (9.6) 

It is trivial to construct two distinct base states (for < v < 2) with energy and 
angular momentum eigenvalues: 

^ o) = n (** - 

a</3 

1 2 



jy = -WV(7V - 1), 



v4°;=n^-%) 2 ^ 

Q</3 

A = -^^-i)- 

The general strategy of the ladder operator approach is to construct multi-anyon 
wavefunctions by acting with step operators on the base states in (10. a) and (10. b). 
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We must, however, respect the statistics and hard-core requirements for the resulting 
wavef unctions. In order to respect the statistics, we use only symmetric combinations 
of the step operators. Consider the symmetric operators 

Q n = X><Ua > (11) 
a=l 

where l,n are non-negative integers such that I + n < N. (These operators form a 
basis in the ring of symmetric polynomials in 2iV variables.) They are step operators 
in energy and angular momentum which respect the statistics properties of the base 
states 

[H, Ci n ] = ulCi n , (12.a) 
[J,C ln ] = (n-l)C ln . (12.6) 

They may, however, produce singular states (states with non-vanishing wavefunctions 
at points of coincidences) which have to be excluded by hand. We must identify which 
particular C\ n produce regular states. These operators can be safely applied to the 
base states. Consider C 0n first. With (8.b), we have 

N 

CWf = (X^f. (13) 

a=l 

Thus they can be safely applied. For Ci m , we have 

TV 

a=l 

N 

= Y,(^ - d a )z™j>f> (14) 

N N 
a—1 a—1 

where we set eB = 2 for simplicity. The seemingly singular first term is in fact regular 
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because 



N N 



Z, 



,m 
a 



— Z 



V 





(15) 



a=l a=l f3^a 



Za ~ Zp 



a<0 



Za ~ Zp 



By a similar proof, one can apply a sequence of operators of the form Co n followed by 
a sequence of operators of the form C\ m to tpP without generating any singularities. 
Moreover, states of the form 



because such a string contains at most derivatives of order Yl)=i n j w hh respect to 



Thus we see that under suitable conditions, the step operators C\ n can be safely 
applied to the base state ^P to generate regular new wavef unctions. A similar 
analysis holds for the other base state ^p ■ Furthermore, closed-form eigenfunctions 
generated by the action of combinations of operators Cn, C\q and Co m only have been 

rgi 

found, and they can be expressed in terms of the Laguerre functions. In particular, 
they do not involve the operators C\ m with m > 1, which however are allowed to act 
on to produce regular wavefunctions. 

One should also note that the step operator approach only generates a subset 

rgi 

of the whole spectrum of wavefunctions. If we naively set v to be zero or one, we 
obtain only a subset of the bosonic and fermionic wavefunctions. Unlike the states 
generated by the step operators, the energies of the missing states show non-linear 
dependence on the statistical parameter v in recent numerical studies.' 61 However, 
this is unimportant for what follows. 




(16) 



a</3 



are obtained from the action of C on on ■?/>} . We can apply a string of operators 
C nimi Cn 2 m2 ■■■ C nimi with Yl)=i n j < 2 ^ to without generating singularities, 
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3. Non-Abelian Chern-Simons Particles 



Recently, there has been much interest in the non-Abelian generalization of 
anyons. Non-Abelian Chern-Simons (NACS) particles carry non-Abelian charges and 
interact with each other through the non-Abelian Chern-Simons term. It has been 
argued that they may have applications in the fractional quantum Hall effect.' 131 Con- 
sider a system of N particles each of which carries a statistical charge corresponding 
to a representation R[ a , a = 1, . . . , TV of a non-Abelian gauge group, which for definite 
we take to be G = SU(2). In the holomorphic gauge, the dynamics of N free SU{2) 
NACS particles is governed by the Hamiltonian' 141 [12] 



N | 



a=l 



h ^— ' 7- — 7a 



d 



Za dz a ' 

where k, a positive integer, is a parameter of the theory and Tg are the SU(2)- 
generators in the representation R[ a . 

The wavefunctions take values in the tensor product of these representations. 

e R k <g) ...®R lN . (18) 

We expand the single-valued wavefunction ip in terms of the conformal blocks Fi e 
R h <g) ■ ■ ■ <S> Ri N (which satisfy V a F; = 0): 

m = Y,^iFi. (19) 

i 

The Hamiltonian acting on the new wavefunction is just the free Hamiltonian. How- 
ever, the complexity of the problem is hidden in the multivaluedness of the wavefunc- 
tions ipi- (In fact, it is more "natural" to work with the multi- valued wavefunctions 
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ipi than the original single- valued wavefunctions ip, partly because in the holomorphic 
gauge the Hamiltonian is not hermitian with respect to the usual inner product. In- 
stead, the inner product is defined in the singular gauge and transformed back to the 
holomorphic gauge by a non-unitary transformation function which has to be taken 
into account in the definition of the inner product.' 121 ) 

From now on, we stick to the singular gauge. Consider N NACS particles in 
the same irreducible representation R[ of SU (2) moving in a uniform external mag- 
netic field B. We introduce operators 

O'oli ^1) b a , ba as in eqn.(8) of section 2 and 
find that the Hamiltonian is again given by eqn.(9). The only difference lies in the 
constraints of the monodromy properties of the wavefunctions. In the case of anyons, 
the wavefunctions have only one component and monodromy leads to acquisition of 
phases, whereas NACS particles have multi-component wavefunctions whose mon- 
odromy properties are given by matrices. 

We define 

a 

Note that J2 a </3 ^a(3i J an d H commute with each other and are thus good quantum 
numbers. (J — J2 a <i3^a/3 is the angular momentum operator in the holomorphic 
gauge.) We will discuss the diagonalization of J2 a <f3^af3 later. For the time being, 
let us assume this has been done and let ipj G <S> Ri 2 <£>...<£> Ri N be a (position- 
independent) eigenvector of J2 a </3 with eigenvalue Q. In analogy with the anyon 
case, we propose applying the same ladder operator approach with the following base 

[15] 

states which are expressed as path-ordered line integrals: 

i)f\zi, . . . ,z N ) = Pexpf / - 2m a pI)d\og(z a - zp) Hi, (21) 

^ a</3 ' 

$n (zi, ...,z N ) = Pexp( / ^2(2n a/3 I - Q a p)d\og(z a - zp) Wj, (22) 

where T is a path in the iV-dimensional complex space with one end point fixed 
and the other being ( = (z\, . . . , zn). The m a p {n a p) depend on ipj and are the 
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maximal (minimal) integers which make the wavefunctions non-singular at the points 
of coincidences. This is analogous to the requirement < v < 2 in the anyon 
case. Modulo the terms involving the identity matrix, the first integrand is just 
the flat Knizhnik-Zamolodchikov connection 1161 whereas the second is related to its 
antiholomorphic analogue. One can easily check that these base states have the 
desirable monodromy properties. From (5), (6) and (7. a), we have 



H$f> = 0, 



E°j = -Nu, 

a</3 

= (fi-2^m Q/3 )^ 0) , 

a<(3 

a</3 



(23) 



and 



N D 

~m Za 2s 

a=l f3^a 



(2n a/3 I - VL aj3 ) w ) 
yjj 



a<(3 
eB 
m 



(0) 

n 



a</3 



(0) 
II 
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Q=l /3^Q 

= [-2 ^n^ + Q]^ 1 ?, 

Q</3 



(24) 



where we have used the relation 

[J2^a0,^s] = 0, (25) 

and the fact that ipj is an eigenstate of the operator YlaKd^aB- Mathematically, 
these commutator relations are just consequences of the integrability condition (in- 
finitesimal pure braid relations) satisfied by the connection. Physically, they follow 
from the fact that Q is related to the angular momentum J which is invariant upon 
monodromy. 

Now that we have found the analogous base states, we will apply the ladder op- 
erators to them to generate new states. As before, the new states have to respect the 
statistics. (The NACS particles in the same irreducible representation are regarded as 
indistinguishable. f 7] Thus, we may only use symmetric combinations of step opera- 
tors. Also, we have to check that the wavefunctions produced are regular at the point 
of coincidence. There is, however, one crucial difference between the cases of anyons 
and NACS particles. Even with symmetric step operators, there is no guarantee that 
the monodromy properties of the wavefunctions are preserved. Any combination of 
step operators which do not preserve the monodromy properties of the wavefunctions 
are to be rejected. 

First of all, let us consider Cq u . As before, we get 

Co4 M£«k 0) . (26) 
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This shows that Cq h can be safely applied to ijjf 1 ^ without changing its monodromy 
properties or producing singularities. Next we consider C\ m . 

N 

a=l 
N 

= - E WcdT + - ™%r^T (27) 

a=l a=l 



a</3 H a=l 



For m — 0, 



N 
a=l 

which clearly preserves the monodromy property of . When m — 1, we have 

TV 

di^f = - - 2 ^/)^ 0) + 5>«^ - 

a</3 a=l 
TV 

= [-n + 2 ^ m Q/3 + E(*a*a - 1)]^ 0) - 



(29) 



This shows that Cn can be safely applied to the base state. We can say more: strings 
made up of combinations of the operators C$n ( n= l)2,...), C\q and C\\ act on the base 
state to generate physical states. On the other hand, the operators C\ m with m > 1 
and C nm with n > 1 generally change the monodromy property of base state. There 
is no obvious way of constructing an admissible combination of operators involving 
them which would preserve the monodromy property of the base state. We therefore 
reject them as being unphysical and restrict the admissible set of operators to those 
generated by Co n ,Cio and C\\. 

The crucial reason why the argument for C\ m (with m > 1) and C nm (with 
n > 1) as physical operators for anyons do not carry over to the case of non-Abelian 
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C-S particles is that the the various monodromy matrices do not commute. In other 
words, 

[n aP ,n lS ]^o. (30) 



As in the anyon case, there are again missing states in the spectrum. However, our 
wavefunctions do cover the entire lowest Landau level as they involve the operators 
C 0n only. 

We now consider the construction of closed-form eigenfunctions. In the case of 
anyons, for P(z\, . . . , zn, z\, . . . , zn) Yl a <p( z a ~~ z pY to be an eigenfunction of H, it 
follows that the function P has to satisfy a modified differential equation. 

£ Ll 9a + f^W - 1 E (hz*L\ vP = (s - 4) p (3D 

An ansatz has been made to construct closed-form eigenfunctions. All the solutions 
constructed can be expressed in terms of the Laguerre Polynomials. They are gen- 
erated by strings of operators Co n ,Cio andCn only. For NACS particles, we get a 
similar equation for P, but with v replaced by Q a p — 2m a p. Nevertheless, since these 
operators are chosen to preserve the monodromy properties of the states. There is 
every reason to believe that the construction of closed-form solutions will go through 
with ^iV(iV — \)v replaced by VL — 2 X] Q </3 m a/3- 

Finally, we come to the diagonalization of J2 a </3 Consider the identity 

N 

^(Tf + T 2 a + . . . + Tfr)(T? + T 2 a + . . . + T%) = T a T a + 2 Y1 T * T P- ( 32 ) 

a a=l a</3 

For SU(2), the left-hand side gives the Casimir operator J(J + 1) of the "spin" of 
the composite made up of the N particles, and the first term on the right-hand 
side gives the sum of the Casimir operators Yla=i Jai^a + 1) of the "spins" for the 
individual particles. (Here we abuse the word "spin" for the internal SU (2) symmetry 
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group. The physical spin (which is a scalar in 2+1 d) of a NACS particle in the j 
representation is given by | J(J + 1). Thus, for SU(2), we have 

1 N 

n = J2^P = ^[J(J + 1) - J -( J - + !)]• ( 33 ) 

a</3 a—1 

We just decompose the composite state into irreducible representations and 
J2 a< /3^a/3 would be diagonal in that basis. Actually, we can do better than that. 
It is easy to check that the operator Tf + T| + . . . + Tfo commutes with H, J and 
J2 a </3 Thus they can be simultaneously diagonalised. 

4. Second virial coefficient and the large k limit 

In this section, we compute the second virial coefficients for some simple systems 
of NACS particles. To do so, we need to know all the two-particle states only. 

First of all, consider two identical NACS particles in the j — \ representation of 
SU(2). From the addition rule for angular momenta, we find that the resulting states 
consist of a triplet with f2 = and a singlet with f2 = — Jr. For N — 2, Q plays the 
role of the anyon phase, v. Let us recall the formula derived by Arovas et al. for the 
second virial coefficient of anyons, 

B(u = 2j + 5,T) = 4(-\ + \S\- 1 -5 2 ), (34) 

where \5\ < 2. Note that it has a cusp at Bose values v = 2j. By taking the average 
over the four two-body states, the second virial coefficient of the NACS particles is 
given by 

B0 = 5,r) = Aj[-i + l-^]. (35) 

For two particles with j = 1, the resulting states have "spins" 2, 1, and (with 
Q — |, — |, and — | and degeneracies 5, 3 and 1 respectively). We remark that all 
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these states are bosonic if k — 1. When k = 2, the singlet is a bosonic state whereas 
others are fermionic. For k > 4, all the states are anyonic with |z/| < 1. For k > 1, 
we have 

B = l,r) = A3.[-i + |-^]. (36) 
Now we come to the large /c limit. For two particles belonging to a representation 

•2 

j with limj^oo j- — a < 1, we approximate the sum over all the resulting "spins" 

fc^oo 

r < 2j by an integral. For example, the \5\ term is given by 

2i 



r(r 2 -2f) (37) 



r=0 

2 



= J_ 

k ' 

where in the second line we approximate the sum by an integral and divide it into 
two parts (which happen to be equal) according to the sign of r 2 — 2j 2 . The 5 2 term 
can be evaluated in a similar manner. Hence, we get 

B = ^-i + T-^ = ^-I + -Tl- (38) 



If a — > 0, the last term may be discarded and the second virial coefficient of the 
NACS particle in the j representation (with physical spin |j(j + 1)) is the same as 
that of an anyon with half of the physical spin as its statistical parameter. 
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5. Concluding Remarks 

(1) The iV-free-NACS-particle problem can be solved by a similar method. The 
free Hamiltonian in the "anyon" gauge is given by 

N 2 

H = J2 9 a d a . (39) 

a=l 

The subtlety is that our base states become unnormalizable! 91 Let us define r = 
(52a=l \ z a\ 2 )^ an d consider 

HM(r)i)f ] = ijf ] x --^-[d 2 /dr 2 + (l/r)(2JV - 1 + 2Q - 4 ^ m a/3 )]M(r). (40) 

a<(3 

We have eigenfunctions of the form 

M^r) = r-^J^kr), (41) 

(where \x = 2(N — 1 + f2 — 2 2~2a</3 m af3)) with eigenvalues hk 2 /2m for H. 

(2) Let us now consider the construction of Ci m for multi-species non-Abelian 
C-S particles (particles in various irreducible representations). In this case, when we 
construct Co n , we do so for each irreducible representation R and symmetrize over 
particles in this irreducible representation only' 101 Let us call the resulting operator 
Cq u . If we construct and C-q in a similar manner, we find that these operators 
have to be rejected: They do not preserve the monodromy properties of the base 
states because [Q a p,Q^s] ^ 0, and we are no longer summing over all the particles. 
Therefore, for Cio and C\\ we do sum over all the particles in the various irreducible 
representations. 

(3) Note that Coi and Cio represent center of mass excitations. The operator C\q 

[18] 

was also analyzed by Johnson and Canright, while C\\ is directly related to the Lie 
group generator of SU(1, 1). [191 
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(4) We remark that the operators of C\ m and C nm (m, n > 2) do preserve the 
monodromy property of the base state, if ipi in eqns.(20) and (21) is chosen to be a 
simultaneous eigenstate of all Q a 8- The statistics is "Abelianized" in this case. This 
situation occurs, for example, for some models of non-Abelian vortices of finite gauge 
groups such as the quarternion group. 

(5) The same ladder operator approach may well apply to non-abelian vortices of 
finite gauge groups.' 171201 Unfortunately, we generally do not know how to construct 
connections which would produce the desirable monodromy in this case. 

[21] 

(6) In a recent paper, Dasnieres de Veigy and Ourvy derived the equation of 
state of an anyon gas in a strong magnetic field at low temperatures. The idea is that 
at sufficiently low temperatures, excitations to higher Landau levels can be neglected. 
Thus one may consider only the lowest Landau states of the anyons, which are covered 
by the step operators. In fact, apart from the statistical phase factor, the multi- anyon 
states in the lowest Landau level are tensor-product states of the individual anyon 
states. By regularizing the grand partition function with a harmonic potential, the 
equation of state can be obtained. The same decoupling principle should apply to 
NACS particles. For a fixed base state, modulo the statistical term involving Q a p, 
the multi-particle wavefunctions in the lowest Landau level are again tensor products 
of individual particle states. Therefore, in principle, one should be able to derive the 
equation of state of NACS particles in a strong magnetic field at low temperatures. 

We thank John Preskill and Piljin Yi for useful discussions. 
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